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We study 2-field inflation models based on the "large-volume" flux compactification of type 
^vq ' IIB string theory. The role of the inflaton is played by a Kahler modulus r corresponding to a 

4-cycle volume and its axionic partner 8. The freedom associated with the choice of Calabi-Yau 
manifold and the non-perturbative effects defining the potential V(t, 8) and kinetic parameters of 
the moduli brings an unavoidable statistical element to theory prior probabilities within the low 
energy landscape. The further randomness of (t, 8) initial conditions allows for a large ensemble 
of trajectories. Features in the ensemble of histories include "roulette trajectories" , with long- 
s' ' lasting inflations in the direction of the rolling axion, enhanced in number of e-foldings over 
those restricted to lie in the r-trough. Asymptotic flatness of the potential makes possible an 
eternal stochastic self-reproducing inflation. A wide variety of potentials and inflaton trajectories 
agree with the cosmic microwave background and large scale structure data. In particular, the 
observed scalar tilt with weak or no running can be achieved in spite of a nearly critical de Sitter 
deceleration parameter and consequently a low gravity wave power relative to the scalar curvature 
power. 



1 Introduction 

The "top-down" approach to inflation seeks to determine cosmological consequences beginning 
with inflation scenarios motivated by ever-evolving fundamental theory. Most recent attention 
has been given to top-down models that realize inflation with string theory. This involves the 
construction of a stable six-dimensional compactification and a four-dimensional extended de 
Sitter (dS) vacuum which corresponds to the present-day (late-time) universe e.g., the KKLT 
prescription pQ. Given this, there is a time-dependent, transient non-equilibrium inflationary 
flow in four dimensions towards the stable state, possibly involving dynamics in both sectors. 

Currently, attempts to embed inflation in string theory are far from unique, and indeed 
somewhat confused, with many possibilities suggested to engineer inflation, using different axionic 
and moduli fields [2j [3] , branes in warped geometry [4] , D3-D7 models [5] , etc. [7] . These 
pictures are increasingly being considered within a string theory landscape populated locally by 
many scalar fields. 

Different realizations of stringy inflation may not be mutually incompatible, but rather may 
arise in different regions of the landscape, leading to a complex statistical phase space of solutions. 
Indeed inflation driven by one mechanism can turn into inflation driven by another, e.g., [8], 
thereby increasing the probability of inflation over a single mechanism scenario. 

So far all known string inflation models require significant fine-tuning. There are two classes 
that are generally discussed involving moduli. One is where the inflaton is identified with brane 
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inter-distances. Often the effective mass is too large (above the Hubble parameter) to allow 
acceleration for enough e-folds, if at all. To realize slow-roll inflation, the effective inflaton mass 
should be smaller than the Hubble parameter during inflation, m 2 < H 2 . Scalar fields (f> which 
are not minimally but conformally coupled to gravity acquire effective mass terms J2 R(f> 2 which 
prevents slow-roll. An example of this problem is warped brane inflation where the inflaton 
is conformally coupled to the four-dimensional gravity A similar problem also arises in 
supergravity. The case has been constructed that has masses below the Hubble parameter H 
which avoids this ^-problem at the price of severe fine-tuning [3]. Another class is geometrical 
moduli such as Kahler moduli associated with 4-cycle volumes in a compactifying Calabi-Yau 
manifold as in [H [3], which has been recently explored in [9] and which we extend here to 
illustrate the statistical nature of possible inflation histories. 

Different models of inflation predict different spectra for scalar and tensor cosmological fluc- 
tuations. From cosmic microwave background and other large scale structure experiments one 
can hope to reconstruct the underlying theory that gave rise to them, over the albeit limited 
observable range. Introduction of a multiple-field phase space leading to many possible infla- 
tionary trajectories necessarily brings a statistical element prior to the constraints imposed by 
data. That is, a theory of inflation embedded in the landscape will lead to a broad theory "prior" 
probability that will be updated and sharpened into a "posteriori" probability through the action 
of the data, as expressed by the likelihood, which is a conditional probability of the inflationary 
trajectories given the data. All we can hope to reconstruct is not a unique underlying accelera- 
tion history with data-determined error bars, but an ensemble-averaged acceleration history with 
data-plus-theory error bars [TU] , 

The results will obviously be very dependent upon the theory prior. In general all that is 
required of the theory prior is that inflation occurs over enough e-foldings to satisfy our homo- 
geneity and isotropy constraints and that the universe preheats (and that life of some sort forms) 
— and indeed those too are data constraints rather than a priori theory constraints. Everything 
else at this stage is theoretical prejudice. A general approach in which equal a priori theory pri- 
ors for acceleration histories are scanned by Markov Chain Monte Carlo methods which pass the 
derived scalar and tensor power spectra though cosmic microwave background anisotropy data 
and large scale clustering data is described in [lOj . But since many allowed trajectories would 
require highly baroque theories to give rise to them, it is essential to explore priors informed by 
theory, in our case string- motivated priors. 

The old top-down view was that the theory prior would be a delta-function of the correct one 
and only theory. The new view is that the theory prior is a probability distribution on an energy 
landscape whose features are at best only glimpsed, with huge number of potential minima, and 
inflation being the late stage flow in the low energy structure toward these minima. 

In the picture we adopt for this paper, the flow is of collective geometrical coordinates asso- 
ciated with the settling down of the compactification of extra dimensions. The observed inflaton 
would be the last (complex) Kahler modulus to settle down. We shall call this T2. The set- 
tling of other Kahler moduli associated with 4-cycle volumes, Ts,T±, ... and the overall volume 
modulus, Ti, as well as "complex structure" moduli and the dilaton and its axionic partner, 
would have occurred earlier, associated with higher energy dynamics, possibly inflations, that 
became stabilized at their effective minima. The model is illustrated by the cartoon Fig. [T] Wc 
work within the "large volume" moduli stabilization model suggested in [TTJ [T21 [13] in which the 
effective potential has a stable minimum at a large value of the compactified internal volume, 
V ~ 10 5 — 10 20 in string units. An advantage of this model is that the minimum exists for generic 
values of parameters, e.g., of the flux contribution to the superpotential W$. (This is in contrast 
to the related KKLT stabilization scheme in which the tree-level Wq is fine-tuned at ~ 10~ 4 in 
stringy units in order for the V minimum to exist.) 

In this paper, we often express quantities in the relatively small "stringy units" m s oc Mp/VV, 
related to the (reduced) Planck mass 

M P = I/VS^G = 2.4 x 10 18 GeV, (1) 

where G is Newton's constant. 
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Figure 1: Schematic illustration of the ingredients in Kahler moduli inflation. The four-cycles of the 
CY are the Kahler moduli Tj which govern the sizes of different holes in the manifold. We assume T3 
and the overall scale T\ are already stabilized, while the last modulus to stabilize, T2, drives inflation 
while settling down to its minimum. The imaginary parts of Tj have to be left to the imagination. 
The outer 3 + 1 observable dimensions are also not shown. 

In this picture, the theory prior would itself be a Bayesian product of a number of conditional 
probabilities: (1) of manifold configuration defining the moduli; (2) of parameters defining the 
effective potential and the non-canonical kinetic energy of the moduli, given the manifold struc- 
ture; (3) of the initial conditions for the moduli and their field momenta given the potentials. 
The latter will depend upon exactly how the "rain down" from higher energies occurs to populate 
T2 initial conditions. An effective complication occurs because of the so-called eternal inflation 
regime, when the stochastic kicks that the inflaton feels in an e-folding can be as large as the clas- 
sical drift. This T2-model is in fact another example of stringy inflation with self-reproduction. 
(See [2] for another case.) If other higher-energy moduli are frozen out, most inflationary trajec- 
tories would emerge from this quantum domain. However we expect other quantum domains for 
the higher-energy moduli to also feed the T2 initial conditions, so we treat these as arbitrary. 

The Kahler moduli are flat directions at the stringy tree level. The reason this picture works 
is that the leading non-perturbative (instanton) and perturbative (a') corrections introduce only 
an exponentially flat dependence on the Kahler moduli, avoiding the ^-problem. Conlon and 
Quevedo [9] focused on the real part of T2 as the inflaton and showed that slow-roll inflation with 
enough e-foldings was possible. A modification [14] of the model considered inflation in a new V 
direction but with a negative result. 

The fields Ti are complex, Ti = Ti + i9{. In this paper we extend the model of [9] to include 
the axionic direction 9 2 - There is essentially only one trajectory if 9 2 is forced to be fixed at 
its trough, as in |9J. The terrain in the scalar potential K(t2,#2) has hills and valleys in the O2 
direction which results in an ensemble of trajectories depending upon the initial values of T2, #2- 
The field momenta may also be arbitrary but their values quickly relax to attractor values. The 
paper [15] considered inflation only along the 9 direction while the dynamics in the r direction 
were artificially frozen. We find motion in r always accompanies motion in 9. 

In Kahler moduli models, there is an issue of higher order perturbative corrections. Even 
a tiny quadratic term would break the exponential flatness of the inflaton potential and could 
make the 7/-problem reappear. However, the higher order terms which depend on the inflaton 
only through the overall volume V of the Calabi-Yau manifold will not introduce any mass terms 
for the inflaton. Although these corrections may give rise to a mass term for the inflaton, it might 
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have a limited effect on the crucial last sixty e-folds. 

In § [2] we describe the model in the context of type IIB string theory. In § [3] we address 
whether higher (sub-leading) perturbative corrections introduce a dangerous mass term for the 
inflaton. In §0]we discuss the effective potential for the volume, Kahler moduli and axion fields, 
showing with 3 moduli that stabilization of two of them can be sustained even as the inflaton 
T2 evolves. Therefore in §[5] we restrict ourselves to V(t2,&2) with the other moduli stabilized 
at their minima. §[6] explores inflationary trajectories generated with that potential, for various 
choices of potential parameters and initial conditions. In §[7] we investigate the diffusion/drift 
boundary and the possibility of self-reproduction. In § [S] we summarize our results and outline 
issues requiring further consideration, such as the complication in power spectra computation 
that follows from the T2,#2 freedom. 



2 The Type IIB String Theory Model 

Our inflationary model is based on the "large- volume" moduli stabilization mechanism of QTJ [T^] 
[13] • This mechanism relies upon the fixing of the Kahler moduli in IIB flux compactifications 
on Calabi-Yau (CY) manifolds M by non-perturbative as well as perturbative effects. As argued 
in [TTJ 021 H3] , a minimum of the moduli potential in the effective 4d theory exists for a large 
class of models. The only restriction is that there should be more complex structure moduli in 
the compactification than Kahler moduli, i.e. h 1 ' 2 > h 1 ' 1 > 1, where /i 1 ' 2 , ft 1 ' 1 are the Hodge 
numbers of the CY. (The number of complex structure moduli is h ' and the number of Kahler 
moduli is h 1 ' 1 . Other Hodge numbers are fixed for a CY threefold.) The "large-volume" moduli 
stabilization mechanism is an alternative to the KKLT one, although it shares some features with 
KKLT. The purpose of this section is to briefly explain the model of [TT J \V2 \ IT3]. 

An effective 4d Jsf = 1 supergravity is completely specified by a Kahler potential, superpo- 
tential and gauge kinetic function. In the scalar field sector of the theory the action is 



SV=i = / d A x\ 

where 



M 
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- ic^d^d^ - v(M) 



(2) 



V(cj), <f>) = e K/M P ( Kfi DiWDjW - -ttjWW ) + D-terms. (3) 
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Here K. and W are the Kahler potential and the superpotential respectively, Mp is the reduced 
Planck mass eq.(UJ), and 4> l represent all scalar moduli. (We closely follow the notations of [13] 
and keep Mp and other numerical factors explicit.) 
The a' 3 -corrected Kahler potential pS] is 




- ln(5 + S) - In ( -i j O A Q\ . (4) 

Here V is the volume of the CY manifold Ad in units of the string length l s = 2ir^/a' ', 
V = V s ^ and we set a' = 1. The second term in the logarithm represents the a'-corrections with 
£ = ~ ^2(2rr)^ proportional to the Euler characteristic x(M) of the manifold M. S — —iCo + e~^ 
is the IIB axio-dilaton with <f> the dilaton component and Co the Ramond-Ramond 0-form. CI is 
the holomorphic 3-form of M . The superpotential depends explicitly upon the Kahler moduli Ti 
when non-perturbative corrections are included 

W= 9 -l^l{wo + TA l e-^}, Wb = ~ / r T ',A<>. 



_ — 1 ' * \j 1 / -11-%^ , "u — 79 

4vr \ f-J / lj 



Here, Wo is the tree level flux-induced superpotential which is related to the IIB flux 3-form 
Gz = F 3 ~iSH 3 as shown. The exponential terms Aie~ aiTi are from non-perturbative (instanton) 
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effects. (For simplicity, we ignore higher instanton corrections. This should be valid as long as 
we restrict ourselves to 3> 1, which we do.) The Kahler moduli are complex, 

Ti = n + i9i , (6) 

with Ti the 4-cycle volume and 9i its axionic partner, arising from the Ramond-Ramond 4-form 
C4 . The At encode threshold corrections. In general they are functions of the complex structure 
moduli and are independent of the Kahler moduli. This follows from the requirement that W 
is a holomorphic function of complex scalar helds and therefore can depend on Ti only via the 
combination Tj = Ti + i6i. On the other hand, W should respect the axion shift symmetry 
8i — > Oi + — anc j th us cannot be a polynomial function of 7*. (See [17] for discussion.) 

The critical parameters in the potential are constants which depend upon the specific 
nature of the dominant non-perturbative mechanism. For example, a$ = 2zn for Euclidean D3- 

brane instantons and dj = for the gaugino condensate on the D7 brane world-volume. We 
vary them freely in our exploration of trajectories in different potentials. 

It is known that both the dilaton and the complex structure moduli can be stabilized in a model 
with a tree level superpotential Wq induced by generic supersymmetric fluxes (see e.g. |18j ) and 
the lowest-order (i.e. £ = 0) Kahler potential, whereas the Kahler moduli are left undetermined 
in this procedure (hence are "no scale" models). Including both leading perturbative and non- 
perturbative corrections and integrating out the dilaton and the complex structure moduli, one 
obtains a potential for the Kahler moduli which in general has two types of minima. The first type 
is the KKLT minima [1] which requires significant fine tuning of Wq (~ 10 ) for their existence. 
As pointed out in [13) , the KKLT approach has a few shortcoming, among which are the limited 
range of validity of the KKLT effective action (due to a' corrections) and the fact that either 
the dilaton or some of the complex structure moduli typically become tachyonic at the minimum 
for the Kahler modulus. (We note, however, that [TH] argued that a consistent KKLT-type 
model with all moduli properly stabilized can be found.) The second type is the "large-volume" 
AdS minima studied in [111 1121 113) . These minima exist in a broad class of models and at 
arbitrary values of parameters. An important characteristic feature of these models is that the 
stabilized volume of the internal manifold is exponentially large, Vmin ~ exp (ar m m), and can be 
O(10 5 — 10 20 ) in string units. (Here T m in is the value of r at its minimum.) The relation between 
the Planck scale and string scale is 



P = 2 m s ' ( 7 

9 s 



where V m in is the volume in string units at the minimum of the potential. Thus these models 
can have m s in the range between the GUT and TeV scale. In these models one can compute 
the spectrum of low-energy particles and soft supersymmetry breaking terms after stabilizing all 
moduli, which makes them especially attractive phenomenologically (see [T31 120] ). 

Conlon and Quevedo [9] studied inflation in these models and showed that there is at least 
one natural inflationary direction in the Kahler moduli space. The non-perturbative corrections 
in the superpotential eq.© depend exponentially on the Kahler moduli Ti, and realize by eq.© 
exponentially flat inflationary potentials, the first time this has arisen from string theory. As 
mentioned in §[U higher (sub- leading) a' and string loop corrections could, in principle, introduce 
a small polynomial dependence on the Tj which would beat exponential flatness at large values 
of the Ti. Although the exact form of these corrections is not known, we assume in this paper 
that they are not important for the values of the Tj during the last stage of inflation (see § [3]) . 

After stabilizing the dilaton and the complex structure moduli we can identify the string 
coupling as g s = e^, so the Kahler potential ((H) takes the simple form 



M2 =-2 In I V + ^ ) +lng s + K cs , (8) 



where K. cs is a constant. Using this formula together with equations |3]), and (fTTj) . one can 
compute the scalar potential. In our subsequent analysis, we shall absorb the constant factor 
e K-cs m t the parameters Wq and A t . 
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The volume of the internal CY manifold M can be expressed in terms of the 2-cycle moduli 
t\ i = 1, ...,n= h 11 : 

V = ^K ijk tHH k , (9) 
where Kijh is the triple intersection form of M. The 4-cycle moduli Tj are related to the t % by 

9 v 1 



V = -n v]k tH k , (10) 



which gives V an implicit dependence on the r^, and thus /C through eq.®. It is known [21] that 
for a CY manifold the matrix gpgp has signature (1, ft 1 ' 1 — 1), with one positive eigenvalue and 
ft, 1 ' 1 — 1 negative eigenvalues. Since = Ti(P) is just a change of variables, the matrix q®.q t . 
also has signature (1, h 1 ' — 1). In the case where each of the 4-cycles has a non- vanishing triple 
intersection only with itself, the matrix ^..q t . is diagonal and its signature is manifest. The 
volume in this case takes a particularly simple form in terms of the n : 

V = a (V 2 - f>f/j . (11) 

Here a and A, are positive constants depending on the particular model0 This formula suggests 
a "Swiss-cheese" picture of a CY, in which n describes the 4-cycle of maximal size and r 2 , . . . , r„ 
the blow-up cycles. The modulus t\ controls the overall scale of the CY and can take an arbitrarily 
large value, whereas t 2 , . . . , r n describe the holes in the CY and cannot be larger than the overall 
size of the manifold. As argued in H3], for generic values of the parameters Wo, Ai, at one 
finds that n 3> t, and V 3> 1 at the minimum of the effective potential. In other words, the sizes 
of the holes are generically much smaller than the overall size of the CY. 

The role of the inflaton in the model of [5] is the last modulus among the Tj, i = 2, ...,n, 
to attain its minimum. As noted by [9] , the simplified form of the volume eq. (jlip is not really 
necessary to have inflation. For our analysis to be correct, it would be enough to consider a model 
with at least one Kahler modulus whose only non-zero triple intersection is with itself, i.e., 

V = a(...-A 2 rf /2 ), (12) 

and which has its own non-perturbative term in the superpotential eq.©. 



3 Perturbative Corrections 

There are several types of perturbative corrections that could modify the classical potential on 
the Kahler moduli space: those related to higher string modes, or a'-corrections, coming from 
the higher derivative terms in both bulk and source (brane) effective actions; and string loop, or 
(^-corrections, coming from closed and open string loop diagrams. 

As we mentioned before, a'-corrections are an important ingredient of the "large volume" 
compactification models of [TTJ [TH H3] • They are necessary for the existence of the large volume 
minimum of the effective potential in the models with Kahler moduli "lifted" by instanton terms 
in the superpotential. The leading a'-corrections to the potential arise from the higher derivative 
terms in the ten dimensional IIB action at the order ~ a' 3 , 

5 /JB = --V / d w x v/^P)e- 2 '[iZ + 4(^) 2 + a' 3 -^ J - a' 3 Q + ...}, (13) 

ZK 10 J 6-1 4 

where J ~ (R MN pq) a and Q is a generalization of the six-dimensional Euler integrand, 

f d e x^gQ = X . (14) 
Jm 

1 For example, the two-Kahler model with the orientifold of Pf 1>l>li8i g] studied in (22 [H [13] has a = 1/9^2, Ai = 1, 
and A2 = — 1. 
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Performing a compactification of (fT5|) on a CY threefold, one finds a'-corrections to the metric 
on the Kahler moduli space, which can be described by the £-term in the Kahler potential ([4]) 
(see j23[[T6]). We will see later that this correction introduces a positive term ~ £W§ /V 3 into 
the potential. As discussed in [13], further higher derivative bulk corrections at 0(a' 4 ) and above 
are sub-leading to the 1/V 3 term and therefore suppressed. (Note that in the models we are 
dealing with, there is effectively one more expansion parameter, 1/V, due to the large value of 
the stabilized V.) Also, a'-corrections from the D3/D7 brane actions depend on 4d space-time 
curvatures and, therefore, do not contribute to the potential. String loop corrections to the 
Kahler potential come from the Klein bottle, annulus and Mobius strip diagrams computed in 
for the models compactified on the orientifolds of tori. The Kahler potential including both 
leading a' and loop corrections can be schematically written as [M] 

^ = - 21nV -| + v^ + v^ + -- (15) 

(We have dropped terms depending only on the brane and complex structure moduli.) Here /i 
and / 2 are functions of the moduli whose forms are unknown for a generic CY manifold. If they 
depend upon the inflaton Tj polynomially, a mass term will arise for Tj with the possibility of 
an ry-problem. Further study is needed to decide. Note that, although the exact form of higher 
(sub-leading) corrections is unknown, any correction which introduces dependence on t$ only via 
V will not generate any new mass terms for Tj. 

As well as T-corrections there are possible ^-corrections. Non-perturbative effects modify 
the superpotential by breaking the shift symmetry, making it discrete. As noted we do include 
these. Although leading perturbative terms leave the Kahler potential ^-independent, subleading 
corrections can lead to ^-dependent modifications, which we ignore here. 

In this paper, we assume that the higher corrections, though possibly destroying slow roll at 
large values of the inflaton, are not important during the last stage of inflation. 



4 Effective Potential and Volume Stabilization 

In this Section, we sketch the derivation of the effective field theory potential starting from equa- 
tions (|3l5l8p . We choose T 2 to be the inflaton field and study its dynamics in the 4-dimensional 
effective theory. We first have to ensure that the volume modulus V and other Kahler moduli 
are trapped in their minima and remain constant or almost constant during inflation. For this 
we have to focus on the effective potential of all relevant fields. 

Given the Kahler potential and the superpotential, it is straightforward but tedious to compute 
the scalar potential as a function of the fields T{. To make all computations we modified the 
SuperCosmology Mathematica package [25] which originally was designed for real scalar fields 
to manipulate complex fields. 

The Kahler potential © gives rise to the Kahler metric Kq = with 

_ 3 Q (4V-g + 6a(EL 2 A fc r fc 3/2 )) _9a l X i X jy /W7 



4(2V + 2 (^+EL 2 A^ 2 )V3 ' « 2(2V + a ' 
K 9A jV T7( a 5 (V + a EL2 A^ 2 )) 1 / 3 3aA,(2V + e + 6aEL 2 A fc r fc 3/2 ) 



2(2V + 2 ' 4(2V + 2 \A7 

This can be inverted to give 

I<n 4(2V + 0(V + aEL 2 ^r fc 3/2 ) 1 / 3 (2V + e + 6ay:L2^ fc 3/2 ) i? _ 8(2V + Qr^- 

3a 4 / 3 (4V-£) ' 4V-£ 

3 8(2V + Qg(g + EL 2 A fc r fc 3/2 )V 3 a 4(2V + gVg(4V - g + 6a EL 2 ^J 2 ) 
4V - C ' 3a(4V - C)A 2 



(17) 
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This is the full expression for an arbitrary number of Kahler moduli T{. The entries of the 
metric contain terms of different orders in the inverse volume. If we were to keep only the 
lowest order terms > 0(pr), the shape of the trajectories we determine in the following sections 
and our conclusions would remain practically unchanged. That is, we are working with higher 
precision than necessary. Note that the kinetic terms for r and 9 are identical, appearing as 
K 22 (8t8t + 8889) in the Lagrangian. 
The resulting potential is 

K n) (4V-£)(2V + £) 2 + 2_, (4V-0(2V + 2 3aA 2 (2V + £) 1 ' 



i=1 

32e- 2a * T *a l A 2 l n(l + a l n) 8W A ie ^ a ^ cos^A) ( 3£ 



E 

i,j = 2 
i<3 



(4V-f)(2V + (4V-0(2V + f) V(2V + 

(4V-0(2V + gV P2(2V + 0(0,7, + a jTj 



4a, t, 



Kpiift = ^ , (19) 



+2(najTiTj) + 24£] + Kpiift ■ 

We have to add here the uplift term Vupiift to get a Minkowski or tiny dS minimum. Uplifting is 
not just a feature needed in string theory models. For example, uplifting is done in QFT to tune 
the constant part of the scalar field potential to zero. At least in string theory there are tools for 
uplifting, whereas in QFT it is a pure tuning (see, e.g., [2 [26] and references therein). We will 
adopt the form 

A 
v 2 

with (3 to be adjusted. 

We now discuss the stabilization of all moduli Tj plus the volume modulus. For this we 
have to find the global minimum of the potential cq. (|18[) . which we do numerically However, it 
is instructive to give analytic estimations. Following [121 113j . we study an asymptotic form of 
eq. ([T5|) in the region where both V ~ exp (djT,-), i = 2, n, and V>1. The potential is then a 
series of inverse powers of V. Keeping the terms up to the order O(^) we obtain 

v _ i £ W-,„ + ^ aiAiWme -„ cosM() + ^ + ,20) 

i—2 i—2 

The cross terms for different Tj do not appear in this asymptotic form, as they would be of order 
1/V 5 . Requiring |^ = and cos(eij0j) = — 1 at the minimum of the potential eq. ([20]) . we get 

V ■ ~ W ° X2 pf^— e ^, min „ W ^ f^-—^,^ „ „ Wok R a„r„, mte (21) 

a 2 A 2 o 3 A 3 a n A n v 

where T i min are the values of the moduli at the global minimum. The expression (|20[) has the 
structure 

where the coefficients C\, C2, C3 are functions of Tj and 9i . C\ and C3 are positive but C 2 can 
be of either sign. However, the potential for the volume V has a minimum only if C 2 < 0, which 
is achieved for cos(a^i) < 0; otherwise V(V) would have a runaway character. Also if all n are 
very large so that e~ aiTi — > 0, then C2 — > and V cannot be stabilized. Therefore to keep C 2 
non-zero and negative we have to require that some of the i > 2 Kahler moduli n and their 
axionic partners are trapped in the minimum. For simplicity we assume all but T 2 are already 
trapped in the minimum. 

It is important to recognize that trapping all moduli but one in the minimum cannot be 
achieved with only two Kahler moduli t\ and t 2 , because t\ effectively corresponds to the volume, 
and t 2 is the inflaton which is to be placed out of the minimum. The Fig. [2] shows the potential 
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(a) 



(b) 



T 



Figure 2: (a) The potential surface V{t\,T2) in a two-Kahler model, with the axionic components 
9\ and #2 fixed at their minima, (b) shows the related contour plot of the volume V against ti- 
Although it may be possible to find a local (very shallow) minimum for the volume in this model 
(marked by a star), the generic situation is that both t\ and t<i will be dynamical, and indeed, the 
evolving T2 could force t\ out of a local minimum, thereby destabilizing what may have once been 
stabilized. For this reason, we have focused on models with three or more Kahler moduli, with all 
but one mutually enforcing their respective stabilizations, and in particular that of the volume. This 
large volume multi-Kahler approach to stabilization differs from the KKLT stabilization mechanism. 

as a function of t± and ti for the two-Kahler model. One can see from this plot that a trajectory 
starting from an initial value for T2 larger than a critical value will have runaway behavior in the 
Ti (volume) direction. Thus, as shown by [9], one has to consider a model with three or more 
Kahler moduli. 

By contrast, the "better racetrack" inflationary model based on the KKLT stabilization is 
achieved with just two Kahler moduli [3]. However, in our class of models with three and more 
Kahler moduli we have more flexibility in parameter space in achieving both stabilization and 
inflation. Another aspect of the work in this paper is that a "large volume" analog of the "better 
racetrack" model may arise. 

We have learned that to be fully general we would allow all other moduli including the volume 
V to be dynamical. This will lead to even richer possibilities than those explored here, where we 
only let T2 evolve, and assume that varying it does not alter the values of the other moduli which 
we pin at the global minimum. To demonstrate this is viable, we need to show the contribution 
of Ti to the position of the minimum is negligible. Following [9], we set all Ti, i — 2, ...,n, and 
their axions 6i to their minima and use equations (j20|) and ([21) to obtain the potential for V: 




(23) 



As one can see from eq. ([2Ul) . the contribution of T2 to the potential is maximal (by absolute 
value) when T2 and 82 are at their minimum, and vanishes as T2 — * +00. This gives a simple 
criterion for whether the minimum for the volume V remains stable during the evolution of T2: 
the functional form of the potential for V (|23|) is insensitive to T2 provided [9] : 




(24) 



For a large enough number of Kahler moduli this condition is automatically satisfied for generic 
values of at and Ai. We conclude that with many Kahler moduli the volume does not change 
during the evolution of the inflaton T2 because the other T stay at their minimum and keep the 
volume stable. 




' B 1.2x10 6 1.4x10 6 



(a) T (b) 

Figure 3: One dimensional sections of the uplifted potential for parameter set 1. We perform a 
proper uplift procedure by explicitly introducing the additional field T3 which is responsible for 
stabilizing the volume during inflation. The parameters for T3 are chosen in such a way that the 
stability condition is fulfilled and at the same time we recover the desired value of V m i n ~ 10 6 . The 
minimum in T3 direction is clearly visible at T3 jm j n ~ 49. 



Consider a toy model with three Kahler moduli in which T2 is the inflaton and T3 stays at 
its own minimum to provide an unvarying minimum for V. We choose parameters as in set 1 in 
Table 1 (which will be explained in detail below in Sec. [5]), and also let 03 = 27r/300, A3 = 1/200, 

and A3 = 10. Eq. (f24)) is strongly satisfied, I -^h = 10~ 4 , under this choice of parameters. 



372 / 3/2 

-2 / a 3 

Therefore we can drop the t 2 , ^-dependent terms in the potential (|20| and use it as a function 
of the two fields V and T3 to find their values at the minimum (after setting also #3 = 71-/0,3 to its 
minimum). The minimization procedure should also allow one to adjust the uplift parameter (3 
in a way that the potential vanishes at its global minimum. With our choice of parameters we 
found the minimum numerically to be at T3 = 49 and V = 10 6 with (3 = 8.5 x 10 -6 , as shown in 
Fig. El 

5 Inflaton Potential 

We now take all moduli Tj, i > 3, and the volume V (hence T\) to be fixed at their minima, but 
let T2 vary, since it is our inflaton. For simplicity in the subsequent sections we drop the explicit 
subscript, setting T2 = t + i8. The scalar potential V(t,8) is obtained from eq. tfTSj) with the 
other Kahler moduli stabilized: 



(4V m -0(2V m +£) 2 (4V m -£)(2V m + £) 2 (2V m +0 



D 3 + 32e- 2a2T a 2 ^T(l + a 2 r) D 4 + 8W A 2 e- a ^ cos(a 2 0) / 3£ \ J3_ 

(4V m - (2 V m + (4V m - (2V m +0 V (2V m + £) ^ ) V? n ' 

Here the terms Di,...,Di contain contributions from the stabilized Kahler moduli other than 
the inflaton. We dropped cross terms between r and other n, i > 3, since these are suppressed 
by inverse powers of V m . We can trust eg. (|25|1 only up to the order 1/V 3 ; at higher orders in 
1/V, higher perturbative corrections to the Kahler potential eq. (fT5"]) start contributing. Explicitly 
expanding to order 1 /V 3 yields the simpler expression 

3a\ 2 V m V£ 
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Parameter 


Wn 




A 2 


A 2 


Oi 


S 


tfs 


V 


Aw/M„ 

— r / P 


Parameter set 1 


300 


27T/3 


0.1 


1 


l/9y/2 

1 v 


0.5 


1/10 

1 


10 6 


2 x 10~ 3 


Parameter set 2 


6 x 10 4 


27r/30 


0.1 


1 


l/9>/2 


0.5 


1/10 


10 8 


1 x 10~ 3 


Parameter set 3 


4 x 10 5 


tt/100 


1 


1 


l/9>/2 


0.5 


1/10 


10 9 


1.4 x 10~ 3 


Parameter set 4 


200 


7T 


0.1 


1 


l/9>/2 


0.5 


1/10 


10 6 


1.5 x 10~ 3 


Parameter set 5 


100 


2vr/3 


0.1 


1 


l/9\/2 


0.5 


1/10 


10 6 


1.9 x 10" 3 


Parameter set 6 


75 


2vr/6 


1 


1 


1/9^2 


0.5 


1/10 


10 8 


4 x 10" 4 



Table 1: Sample parameter sets for the T 2 inflation model. Sets 1 to 4 define models in which the 
primordial power spectra of perturbations are approximately compatible with observations, whereas 
sets 5 and 6 do not. There are approximate scaling relations which map one set of parameters without 
changing the power spectrum (see explanation in the text regarding Wo)- The value of Wq in sets 2 
and 3 is so high that corrections to the potential seem likely to appear. We note however that these 
large values can be reduced using the scaling transformation, but the undesirable cost is that t can 
drop below the string scale or a 2 can become too large. The last column shows approximate values 
for the variation of the canonically- normalized inflaton eq. (|29p over the observable e-fold range, 
appropriately small cf. the Planck scale. 



where 

~ 4V r 3„ + V m + VI {Z1) 

is a constant term, since V and Tj, i > 3 are all stabilized at the minimum, and Z?2, -D4 depend 
only on these rj. 

The potential eq. (j26|) has seven parameters Wo, 02, A2, A2, a, £ , and g s whose meaning was 
explained in § [2j We have investigated the shape of the potential V(t,9) for a range of these 
parameters. Wo, 02, A2 control the low energy phenomenology of this model (see [20]) and are 
the ones we concentrate on here for our inflation application. We shall not deal with particle 
phenomenology aspects in this paper. Some choices of parameters Wo, 02, A2 seem to be more 
natural (see [2Q])- To illustrate the range of potentials, we have chosen the six sets of parameters 
given in Table 1. There is some debate on what are likely values of Wo in string theory. We 
chose a range from intermediate to large. Since there are scaling relations among parameters, 
we can relate the specific ones we have chosen to others. An estimate of the magnitude of Wo 
comes from a relation of the flux 3-forms F3 and H3 which appear in the definition of Wo (eq(5]) 
to the Euler characteristic x °f the F-theory 4-fold, which is \ ~ J F3 A H3 from the tadpole 
cancellation condition. This suggests an approximate upper bound Wo ~ \fx [II]- For typical 
values of x ~ 10 3 , we would have Wo ~ 10 — 100. There are examples of manifolds with \ as large 
as 10 6 , which would result in Wo ~ 10 3 . Further, the bound itself can be evaded by F3 A H3 = 
terms. However, we do not wish to push Wo too high so that we can avoid the effects of higher 
perturbative corrections. We can use the scaling property for the parameters to move the value 
of Wo into a comfortable range. This should be taken into account while examining the table. 

The parameter sets in the table can be divided into two classes: Trajectories in sets 1 ... 4 
produce a spectrum of scalar perturbations that is comparable to the experimentally observed 
one (good parameter sets), whereas trajectories in parameter set 5 and 6 produce spectra whose 
normalization is in disagreement with observations (bad parameter sets). Sets 3 and 4 were 
chosen to large values of Wo to illustrate how things change with this parameter, but we are wary 
that with such large fluxes, other effects may come into play for determining the potential over 
those considered here. A typical potential surface V(r, 9) is shown in Fig. [4] with the isocontours 
of V(t, 9) superimposed. 

The hypersurface V(t,0) has a rich structure. It is periodic in 9 with period 27r/ci2, as seen 
in Fig 2J Along 9 = ^j^ 1 - 1 , where I is integer, the profile of the potential in the r direction is 
that considered in [3]. It has a minimum at some r — r min and gradually saturates towards a 
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(a) (b) 

Figure 4: (a) The T2-potential surface V(r, 9) for the parameter set 1 of Table 1. Equipotential 
contour lines are superposed. Both r and 9 are multiplied by the characteristic scale 02- The 
surfaces for the other parameter sets in the Table are very similar when ensealed, even sets 5 and 6 
with high Wq. In all cases, we manually uplifted the potential to have zero cosmological constant at 
the minimum. The periodicity in the axionic direction and the constancy at large a<2j are manifest. 
If instead of r, we used the canonically- normalized field eq, (|29p which is amplified by VV, the 
undulating nature in 9 at large <p becomes more evident. However, the canonically-normalized 
inflaton is trajectory dependent and not a global function, (b) V(r) for 9 at cos(a2#) = ±1 shows 
how the flow for the positive value from large r would be inward, but the flow would be outward 
for the negative value (and be unstable to 9 perturbations). The dashed line is V^o, the r — ► 00 
asymptote. 
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constant value at large r, V(t) — > I4o(l — Ca2Te~ a2T ), where C is a constant!! Along = 
V(r, 0) falls gradually from a maximum at small r towards the same constant value, V(t) — > 
Voc (1 + Ca2Te~ a2T ). Trajectories beginning at the maximum run away towards large r. Fig. [8]^b) 
shows these two one-dimensional sections of the potential. For all other values of 9 the potential 
interpolates between these two profiles. Thus, at large r, the V surface is almost flat but slightly 
rippled. At small r, the potential in the axion direction is highly peaked. Around the maximum 
of the potential it is locally reminiscent of the "natural inflation" potential involving a pseudo 
Goldstone boson [57] (except that 9 and r must be simultaneously considered), as well as the 
racetrack inflation potential [H [3] . 

There are two approximate scaling symmetries in this model (in the asymptotics of (|26[) for 
V>?), similar to those in [2J, 

ai^u^ai, (Ai,W ,0^u 3 / 2 (Ai,Wo,0 , (nA) -> ufaA) , (A. i; "P s , e) — > (Xi,V s , e) ; 
a^a l7 (A t , Wo) -^v-^ 2 {A u W ) , (n, 9$ (r is 4 ) , V s v 3 V s , (X t , £, e) -> (A,, £, e) . (28) 

which can be used to generate families of models, trading for example large values of Wq for small 
values of r. For instance, applying the u scaling to parameter set 1, we can push the value of 
Wo down to Wo ~ 2, but at the same time pushing r to lie in the range r = 0.1 ... 1.0 during 
inflation, a range which is quite problematic since at such small r higher order string corrections 
would become important. 

More generally, for the supergravity approximation to be valid the parameters have to be ad- 
justed to have T m i n at least a few: r is the four-cycle volume (in string units) and the supergravity 
approximation fails when it is of the string scale. However, even if the SUGRA description in 
terms of the scalar potential is not valid at the minimum, it still can be valid at large t, exactly 
where we wish to realize inflation. The consequence of small T m i n is that the end point of infla- 
tion, i.e. preheating, would have to be described by string theory degrees of freedom. We will 
return to this point in the discussion. 



5.1 The Canonically-normalized Inflaton 

If we define a canonically-normalized field ip by dip 2 /2 = K22 dr 2 , then 

It is therefore volume-suppressed. For inflation restricted to the T2 direction, we identify ip with 
the inflaton. The field change Aip over the many inflationary e- foldings N = — In a/a en d is given 
in the last column in Table 1 for a typical radial (t) trajectory. It is much less than M p . (Here 
the scale factor at the end of inflation is a en d so N goes in the opposite direction to time.) The 
variations of the inflaton and the Hubble parameter wrt N, 

d(p/M *-.V2e, (30) 



dN 
dlnH/M p 



= e,e=l + q, (31) 



dN 

suggest we must have a deceleration parameter q nearly the de Sitter —1 and the Hubble pa- 
rameter H nearly constant over the bulk of the trajectories. This is shown explicitly in § [5] and 
Figures [5] and [HJ The parameter e(N) is the first "slow-roll parameter" , although it only needs 
to be below unity for inflation. With e so small, we are in a very slow-roll situation until near 
the end of inflation when it rapidly rises from approximately zero to unity and beyond. 



2 This type of potential is similar to that derived from the Starobinsky model of inflation |28] with a (i? — jr^iR 2 
Lagrangian via a conformal transformation, V((f>) = 6TvMpM 2 [l — exp( — = )] 2 . 
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Equation (|30[) connects the change in the inflaton field to the tensor to scalar ratio r = Vt/V s > 
since to a good approximation r w 16e. Since we find e<l, we get very small r. The following 
relation 29, 30 gives a lower limit on the field variation in order to make tensor modes detectable 

— « 0.46 J — . 32 
Mpi V0.7 v ; 

We are not close to this bound. If tp/M p is restricted to be < 1 in stringy inflation models, getting 
observable gravity wave signals is not easy. (A possible way out is to have many fields driving 
inflation in the spirit of assisted inflation 

When the trajectory is not in the r direction, the field identified with the canonically- 
normalizcd inflaton becomes trajectory-dependent as we describe in the next section and there 
is no global transformation. That is why all of our potential contour plots have focused on the 
Kahler modulus and its axion rather than on the inflaton. 



6 Inflationary Trajectories 

6.1 The Inflaton Equation of Motion 

We consider a flat FRW universe with scalar factor a(t) and real fields (r, 9). To find trajectories, 
we derive their equations of motion in the Hamiltonian form starting from the four dimensional 
Lagrangian (see [25] and references therein) 

C = V^9 f R + Gijftft '- V) , (33) 



1 ' 

with canonical momentum P, = -J^ = 2a 3 Gj,-0 i , where we used \J—g = a 3 and (pi, i = 1,2 stands 



for (t, 9). (The usual field momentum is Pi /a 3 .) Here the non-canonical kinetic term is 

G ij =K 2 ,S ij , K Q = J^j. (34) 

The Hamiltonian is 

H = P i 4> i -C= -K(>" J I',1', + V , (35) 
where G y = GJ^ 1 . The equations of motion follow from ! = §pp-, Pi — — which reduce to 

1 i r)f kl dV 1 

^ - 2^ ' Pi - - a3 W>" = aH > H = ~^ GijPiP i • (36) 



We also use the constraint equation 

MpH 2 = i (-L&iPiPj+v) =M^H 2 e/3 + V/3 (37) 

to monitor the accuracy of the numerical integration routine. Note that the deceleration, eg. pip , 
is related to the non-canonical kinetic term in the implicit manner indicated. 

The inflaton is defined to be the field combination along the classical (unperturbed) trajectory. 
Isocurvature (isocon) degrees of freedom are those perpendicular to the classical trajectory. The 
kinetic metric in Mp units is dijj 2 = GijdcjfdcfP /Mp. If 9 is fixed, then ip is related to <p introduced 
in § 15.11 by ip = ip/(^/2Mp). More generally, the inflaton between the initial condition <j> l (No) 
and the value <f> l (Ni) is the distance along the path, J dip. For our case with diagonal Gy, it is 



4, 



J v r K^dT(l + d9 2 /dT 2 ) 1/2 . (38) 
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6.2 Stochastic Fluctuations and CMB and LSS Constraints 



The classical trajectory is perturbed by zero point fluctuations in all fields present, but only 
degrees of freedom with small mass will be relevant, hence in r and 9, which in turn influence the 
scalar metric fluctuations encoded in In a = —N, and in the gravitational wave degrees of freedom. 
Structure formation depends upon the fluctuations in the scalar 3-curvature, which are related 
to those in lna measured on uniform Hubble surfaces [32] [33], S^R = A(k / a) 2 5 \n a\ h ■ The 
usual result for single-field inflation is motivated by stochastic inflation considerations: the zero 
point fluctuations in the inflaton at "horizon crossing" when the three-dimensional wavenumber 
k as Ha are Sip = [H/ (2ttMp)]/ y/2. Equation (f3"0|) gives the mapping, (Hna|# = Sip/^Je along 
the inflaton direction. Thus the scalar power spectrum is 

V 8 = k 3 /(2n 2 ) < |<51na| H (fc)| 2 >= [H/(2ttM p )} 2 /(2e) e 2u ' , (39) 

where u s encodes small corrections to this simple stochastic inflation Hawking temperature for- 
mula. 

The graviton zero point oscillations are, like those in a massless scalar field, proportional to 
the Hawking temperature at k = Ha, 

V t (k) = 16 [H/(2nM P )} 2 /2 e 2ut . (40) 

where u t also encodes small corrections. The ratio r(k) is therefore as 16e. 

The scalar spectral index is given by n s — 1 = dlnVs/dlnk. At lowest order, and for small e, 
it is n s — 1 = — 2e — e'/e, where e' = de/dlna. To have e nearly zero as our trajectories do, and 
yet have n s differing at the 2-sigma level from unity as the current cosmic microwave background 
(CMB) and large scale clustering data indicate, imposes a constraint on e'/e which might seem 
to require a fine-tuning of the potential. 

The current best estimate of n s in flat universe models characterized by six parameters, in 
which gravity waves are ignored, is 0.96 ± 0.017 with CMB only, and 0.958 ± 0.015 with CMB 
and large scale structure (LSS) clustering data [33]. The errors are Bayesian 1-sigma ones. For 
future reference, we note that the CMB data prefer a running of the scalar index at about the 
2-sigma level, dn s /d\nk = — 0.047 ± 0.021 at a pivot point k p = 0.05 Mpc -1 . However, the 6 
parameter case with no running is a very good fit, except in the low £ regime. The preferred 
amplitude of the scalar perturbations is V s {k p ) = pl^ 1 ^] x 10 -10 . It is interesting to note that 
this number has been stable for a long time: the estimate from the COBE DMR experiment 
without the addition of any LSS or smaller scale CMB experiments was [21 ± 3] x 10 -10 when 
extrapolated with this n s = 0.96 slope, and only slightly higher with no tilt [33] . 

The current constraint on the gravity wave contribution is Vt/V s < 0.6 at the 95% confidence 
limit with CMB data (with the powers evaluated at the pivot point 0.002 Mpc -1 ). When LSS data 
is added to the CMB, this drops to an upper limit of 0.28 but requires a single slope connection 
of the low k regime in which the tensors can contribute to the CMB signal and high k where the 
amplitude of LSS fluctuations is set. Relaxing this allows for higher values [10] . 

With the CMB-determined V s estimate, we have [H/{2wM P )] as 6.5 x lO -5 ^- If the ac- 
celeration e were uniform over the observable range and gave rise to this n s , we would have 
e = 0.02, [H/(2nMp)] as 10 -5 and r ~ 0.3. But our trajectories have e nearly zero and H 
almost flat, so to get the observed n s the observable range would have to be well into the braking 
period towards preheating: i.e., we would need e'/e as 0.04 over the CMB+LSS window, which 
seems like it would require a very finely-tuned potential. Rather remarkably, the first cases we 
tried gave values near this at the relevant number of e-foldings before the end of inflation. This 
point was also made by Conlon and Quevedo [9] for pure r trajectories. We also find there 
is room for modest running of the scalar index over the observable window. (Note that the 
tensor slope n t = dlnVs/dlnk is as — 2e, hence very small.) A consequence of the small e is 
[H/(2%M P )] - 10 -10 and V 1 ^ - 10 14 Gev. 

We caution that the single field model we used to estimate V s includes fluctuations only along 
the trajectory. There will also be fluctuations perpendicular to it, and fluctuations in the parallel 
and perpendicular directions can influence each other. The latter are isocurvature degrees of 
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freedom. Although these will leave the Vt formula unaffected, we expect modifications in the V s 
formula, a point we return to after making an inventory of the sorts of trajectories that will arise. 

6.3 Trajectories with General Kahler modulus and Axion Initial Con- 
ditions 

6.3.1 The r-valley Attractor 

We first restrict ourselves to trajectories in r to connect with the Conlon and Quevedo [S] treat- 
ment. The stable flow is in the 9 — ^ trough. As can be seen in Fig. [5^a), enough e-foldings 
for successful inflation are possible provided one starts at large enough r. The dashed line in 
Fig. [SJa) shows e is very small for this case, of order e w 1CP 10 . For the parameters we have 
considered, no effective inflation is possible if we start inward of T m i n rather than outward. 

Another class of r trajectories are those along the "ridge" where 9 = 2zl gives a positive 
contribution to the potential. These are unstable to small displacements in the axion direction. 

The r-trough trajectories serve as late-time attractors for initial conditions that begin with 
9 out of the trough. The very flat profile of the potential at large r allows for a regime of self- 
reproducing inflation (§EJ. Trajectories which originate in the self- reproducing regime invariably 
flow to the r-valley attractor and the observed e-folds would be just those of the Conlon and 
Quevedo sort. 

6.3.2 The Variety of "Roulette" t-6 Trajectories 

When we allow the initial values of r and 9 to be populated with an equal a priori probability prior, 
given a set of parameters defining the V(t, 9) surface, we encounter a wide range of inflationary 
trajectories. Examples of the variety of behaviours for the parameter sets given in Table 1 are 
shown for r(ln a), #(ln a) in Figures [5] and [6l and for iJ(lna), e(lna) in Figures [8] and [9] Some of 
the trajectories are predominantly 9 ones before settling into a r-valley, similar to a roulette ball 
rotation before locking into a slot with a specific number, hence the name roulette inflation. 

We began with the momenta of the fields set to zero, but the momenta are very quickly 
attracted to their slow-roll lock-in values, on of order an e-fold. We do not show this settling 
down phase in the trajectories we have plotted. We find a large fraction of trajectories are indeed 
inflating, and have the required > 40 — 50 e-folds of inflation (§ I6.3.4j) to give homogeneity and 
isotropy over our observable Hubble patch. Large enhancements of the number of e-folds over 
r-only inflation can occur because of significant flows in the ^-direction, while r evolves slowly. As 
the figures show, initial values starting far out in r generally roll towards the nearest r-valley and 
then proceed along the r-attractor. Trajectories starting at intermediate r have enough energy 
to pass through the r-valley in the axionic direction, but often turn around before reaching the 
neighbouring ridge, and roll back into the valley. They can move to larger r while in the axion- 
dominated flow. If the initial values arc chosen such that the initial potential energy is just a little 
higher, the inflaton can climb over the next ridge and settle in the adjacent valley, or the next 
one, or the next. Thus there exists bifurcation points which divide the phase space into solutions 
ending up in different valleys. Another feature in Fig. [5] is the existence of areas in which tiny 
changes in initial positions can lead to dramatic changes in the number of e-folds produced. 

There are also cases in which the r-attractor is not reached before the end of inflation. When 
starting from moderate values in r, relatively far away from the valley, the field oscillates in the 
9 direction, albeit producing only a very small number of e-folds N < 1 (see Fig. [5] for examples). 

6.3.3 The Scalar and Tensor Power Spectra and Isocurvature Effects 

We estimate the power spectra V s {k) using the single inflaton approximation, eq. (p?9")) , and Vt{k) 
using eq. (|40[) . which only require H and e for each trajectory. These are shown in Figures [TOl and 
Figures [TTJ As expected from the very small e's we have encountered, to match V s and hence 
H 2 /e to the data would require small H, low energy inflation, and hence very small gravity wave 
power, V t oc [H/M P ] 2 . 
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Figure 5: Contour-plots of the potential including trajectories for several choices of initial values 
(denoted by filled circles) in field space (r, 6). The trajectories are evolved numerically until inflation 
ends at e = 1. The number of e-folds is indicated next to the corresponding trajectory. During the 
last stages of inflation, the field always rolls along one of the valleys towards the minima which 
are located in the centre of the white circles. The maxima are located in the dark spots. Inflation 
in the axion direction can significantly enhance the amount of inflation over that obtained in pure 
r inflation. Trajectories starting at large r roll to the nearest valley, then to the minimum. But 
starting at intermediate r with the axion sufficiently far from its minimum, we find the field can 
cross several #-ridges before settling into a valley. Another manifestation is the run-away character 
for r if the axion is placed close to its maximum, (a) shows the simple pure-r inflation if 9 is set to 
its minimum, as in Conlon and Quevedo, for parameter set 1. (b) shows the complex evolution for 
sample general starting conditions, for parameter set 1, (c), (d) and (e) show the same for sets 2, 3 
and 4, respectively. 
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Figure 6: Sample trajectories for parameter sets 5 (a) and 6 (b), which have scalar power spectra 
amplitudes incompatible with the data. These still look similar to those in Fig. [5j 

For comparison, Figures[lO]and Figures[ll]also show the best-fit scalar power spectra that were 
obtained by the ACBAR collaboration [34] using CMB and LSS clustering data. The fixed power 
law case has n s = 0.96. The spectrum with running has best-fits n s = 0.9f f , dn s /dln k = —0.044. 
Obviously the running model can at best be valid only over a limited number of e-foldings. We 
indicate the observable range of e-foldings by the heavy black line in the power spectrum figures. 
As we show below, the mapping of wavenumber to number of e-foldings is imprecise but because 
the energy of inflation is low, the relevant range for CMB+LSS observables is in the range 
N = — In a e [40, 50]. For definiteness we take the pivot point k p for running to be at N = 45. 

To compute the scalar power spectra we assumed that the generation of fluctuations is driven 
by small perturbations along the trajectories, neglecting the effect of isocurvature fluctuations 
coming from perturbations perpendicular to this direction. We will now discuss why this might 
not necessarily catch the whole picture. In FigJSh 1 , there are a set of trajectories with initial values 
around r s» 21 very close to each other where all fields end up in the second valley from the top. 
Two things are striking about those realizations of inflation: (1) the trajectories flange out when 
rolling down the axion towards larger r's during the first part of their evolution; (2) the number 
of e-folds of inflation varies tremendously from as little as N ss 8 all the way to N 106 when 
going towards larger initial £?'s while keeping r fixed. The scalar spectra we compute for these 
flanging trajectories differ substantially, and make it clear that clear that our simple single field 
algorithm will be unjustified in some regions of the space of initial conditions. Even though the 
transverse quantum jitter is characterized by very small H/Mp with width much smaller than 
the size of the curves in the figure, we recognize that tiny changes during the initial period of 
evolution in certain areas could produce big effects, an area for future investigation. 



6.3.4 Number of e-folds N and the wavenumber of perturbations 

Since our trajectory computations are in terms of In a/a en d, whereas the observables probe k ~ 
Ha in Mpc , we need to connect the two. The CMB+LSS probe k from ~ 10 _4 Mpc to 
~ IMpc -1 , about 10 e-foldings. According to general lore (see e.g. [35]) the number of e-folds 
N as a function of k is given by 

k 

N(k) = 62 - In r + A, (41) 

6.96 x lO-SMpc" 1 

where 6.96 x 10 _s Mpc -1 is the inverse size of the present cosmological horizon and A is defined 
by the physics after inflation: 



10 16 GeV 1 V k 1 
V }I A 4 lr V cnd 3 ln p V4 ' M 

fc " rcn 
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{Number of Efolds:, 84, 109, 2, 28 , 4, 109} 




-3 



Figure 7: Potential with inflationary trajectories for parameter set 4. Shown is a plot of the potential 
in the (r, #)-plane, overlaid with equipotential contours. Note that for r > 5 there are no contours 
because the potential is exponentially flat but it still has the periodic structure in the ^-direction. 
The starting points of the trajectories are indicated by filled circles. 
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Figure 8: Hubble parameter H (In a) (left column) and deceleration parameter e(lna) (right column) 
as a function of the number of e-folds. H(lna) is practically constant during inflation, e during 
inflation is tiny and generally turns up towards the reheating phase rather rapidly, (a) Parameter 
set 1: the coloured trajectories are for general (r, 6) inflation, the dashed trajectory is inflation 
strictly along r (b,c,d) are the same for sets 2, 3, 4, respectively. Note that units along the x-axes 
are In a = —N ~ Ink here and in the plots for the power spectra. There is a very rapid phase in 
which the trajectory settles down to an attractor for the r and field momenta which we do not 
show. 
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Figure 9: Same as Fig. El but for parameter sets 5 (a) and 6 (b) 
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Figure 10: Power spectra for scalar V s and tensor Vt fluctuations, derived assuming that only 
stochastic kicks along the trajectory are relevant for their determination, neglecting the influence 
of "isocurvature" fluctuations transverse to the inflaton trajectory. Both have amplitudes given by 
the instantaneous Hawking temperature. For comparison we present two template spectra: The 
dashed-dotted line shows a simple spectrum with no running of the spectral index n s = 0.96. The 
dotted line is a simple spectrum with running of the spectral index n s = 0.91, dn s /d\nk = —0.044, 
using values obtained by the ACBAR collaboration [33]. Both these spectra have the normalization 
set to V s = 2.1 x 10 -9 at N = 45. (a,b,c,d) are for parameter sets 1, 2, 3 and 4. The dashed line 
corresponds to inflation with 9 =const in the valley. Even though there is significant running for all 
models over large scales, the spectra are mostly featureless in the observable interval which panels 
(e) and (f) zoom into for sets 1 and 4. 
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Figure 11: Power spectra for potentials with parameter sets 5 and 6 that do not give a scalar 
amplitude near to the observed value. 



where p Ten is the energy density at the end of reheating, Vk is the value of the inflaton potential 
at the moment when the mode with the comoving wavenumber k exits the horizon at inflation 
and V cn d is the value at the end of inflation. We will not go into any details about this formula 
or its derivation, but merely motivate some numbers for the individual terms. 

Starting with A, we note that for its first term, H w 3 x 10~ 10 M p i corresponds to V^ 4 ~ 10 13 
GeV. The second term can be neglected in our case, and the last term in the expression for A 
depends on the details of (p)rehcating (which can be perturbative preheating, non-perturbative 
reheating or reheating involving KK-modcs) which we put in the range of 1 < A < 10. Putting it 
all together we find that A 6 [—17, —8]. Therefore the observable range is about 10 e-folds inside 
the interval [35,55], with the exact location of the former depending on the details of reheating. 
For the purpose of our discussion, we take the observable N € [40,50], and indicate this range 
by the black bar in Fig 1101 We note that the required number of e- foldings is significantly lower 
than 65. 



7 Stochastic Regime of Self-Reproduction 

Our potential is flattening exponentially rapidly as r increases. Therefore the velocity of the 
fields which are placed at large enough r will be rather small. This opens the possibility for 
stochastic evolution of the fields due to their quantum fluctuations dominating over classical 
slow-roll [3E1 [32 • We now show that in the model (fT5|) there is region of (r, 9) space where 
the regime of self-reproduction operates. We distinguish this regime from the regime of eternal 
inflation due to bubble nucleations between different string vacua. 

The criterion for self reproduction is actually that scalar perturbations are at least weakly 
nonlinear and that perturbation theory breaks down. The drift in each e-fold of the scalar cp l 
is A(f> 1 = G l ° (Pj/a 3 )/(2H) where Pj/a 3 is the canonically-normalized field momentum. The 
corresponding drift in the normalized inflaton ip is Aip = \fe. The rms diffusion due to stochastic 
kicks is Sip = [H/(2ttMp)]/V2, as given in § 16.21 The rms kick beats the downward drift when 
vV 2 ~ [H/ (2-irMp)}/ V2e exceeds unity, that is the fluctuations become non-perturbative. With 
such a flat nearly de Sitter potential, this is possible. 

We now wish to consider this boundary as a function of 6 as well as r. In general it is 
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Figure 12: The self-reproduction region in which stochastic kicks dominate over the classical down- 
ward drift is indicated by the grey shaded area superposed upon the contour plot of the potential 
and some of its inflationary trajectories for parameter set 1. The most interesting trajectories start 
from outside of the quantum region. Those coming from the stochastic region will be attracted to 
the 9 = r-trough trajectories. 

impossible to bring both of their kinetic terms simultaneously into canonical form. However, 
when we consider the fluctuations of fields, we can revert to the approximation in which the 
Kahler metric stays approximately constant over the time-scales of the fluctuations. The result 
is the following condition for the region where quantum fluctuations dominate over drift: 

v i (v.A 2 v _L(Y±\ 2 

12^ 2 > G TT \ V J ' 12^ 2 > Gg 9 \ V ) ' ( ' 

This result is consistent with starting from (\> % + 3H<j> 1 + T\ k ^P^r + G^V^j = 0, using slow-roll to 

neglect terms of 0{<j) 2 ,<j)) and using the diagonality of G. In Fig. 1 121 we replot potential contours 
with various trajectories for parameter set 1, but with the drift/diffusion boundary now shown. 
Starting trajectories with arbitrary at the boundary should lead to large kicks in 6 as well as 
r, but ultimately a settling into a r-trough well before the observable N range is approached. In 
that case, the 9 complexity of trajectories would not manifest itself. 

Will the universe truly be in a stochastic regime of self-reproduction - the scenario of eternal 
inflation [37] - within such a model? The picture is of the last hole jittering about in size at r 
far from the potential minimum. If the last hole can jitter about far away from equilibrium, then 
the same phenomenon would be expected for other holes that we took to be stabilized before 
the last stages of T-i settle-down. In that case, it is unclear how the initial conditions for T-i 
would be fed, but the most probable source of trajectories would not necessarily be from the Ti 
self-reproduction boundary. In all cases, it is unclear whether further corrections to the potential 
way out there will uplift it to a level in which drift steps exceed diffusion steps. Fortunately, we 
are not dependent on such an asymptotically flat potential for the inflation model explored here 
to work. 

8 Discussion and Summary 

We have investigated inflation in the "large volume" compactification scheme in the type IIB 
string theory model of [TTJ [TSJ [13] . Dynamics in the model are driven by Kahler moduli associated 
with four cycles and their axionic partners, T, = Tj + 9i. We focused on the situation when all 
Kahler moduli/axions and the compact volume are stabilized at their minima except for a last 
modulus/axion pair Ti = (t,0) which we identified with the "observed" early universe inflaton. 
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We showed that Ti dynamics does not perturb the global minimum of the total potential if 
T2 operates on a lower energy scale than that associated with the earlier stabilizations. To do 
this, three or more moduli were needed. We explicitly demonstrated stabilization with a three 
field example, Ti,T2,T%. We also showed volume destabilization can occur even if it is initially 
stabilized if there are only two Kahler moduli, T\,T2- 

The two-dimensional potential V(t, 9) has a rich ridge-trough, hill-valley structure. We de- 
rived and solved the equations of motion of the (r, 9) system appropriate to an expanding universe 
using the Hamiltonian formalism of [25] . The non-canonical nature of the kinetic terms played 
an important role in defining the evolution. 

The ensemble of inflationary trajectories is a rich set. We first studied inflation in the r 
direction, along the valleys of the potential surface, which are late-time attractors for most of 
the general trajectories. For example, trajectories originating in the region of self-reproduction 
wc identified in § 16.21 — where stochastic diffusive kicks can beat classical field drift — would 
invariably finish in the r-troughs. We calculated the number of e-folds N, the time variation of 
the Hubble parameter as a function of In a, -ff(lna), the time variation of the acceleration history 
as encoded in the "first slow roll parameter" e(lna), and the power spectrum for scalar V s and 
tensor Vt fluctuations for many r-trajectories in the different realizations of the potential we 
explored. We can use the tools of single field inflation with confidence to compute the spectra 
if we are only following evolution of the inflaton within the Kahler modulus valley. Our results 
support the conclusions of Conlon and Quevedo [9] who first suggested this model of inflation. 
In particular, trajectories exist in the "prior" ensemble which satisfy the CMB+LSS data. 

Our main objective was to study inflation with general initial conditions for (r, 9). We found 
that trajectories originating in elevated parts of the potential (near its maxima or its ridges) 
with moderate to large initial r values can pass over many valleys and hills in 9 before settling 
into a final r-valley approach to the end of inflation. These roulette wheel trajectories have an 
enhanced numbers of e-foldings relative to pure r trajectories. No fine tuning of parameters 
was needed to get plenty of long-lasting inflation trajectories. In this respect, the Kahler modu- 
lus/axion inflation is superior to "better racetrack" inflation [3|. It is also preferable to have 200 
good trajectories with a single complex field rather than one good trajectory with 200 complex 
fields as occurs in the string version [7] of "assisted inflation" |31) . The power spectra of scalar 
fluctuations for trajectories along a general direction need to be calculated using tools appropri- 
ate to multiple-field inflation. In order to estimate the amplitude we used expression (|39[) for 
single field stochastic inflation. We found indications of weak running in the spectral index over 
the observable range for some trajectories, and no running for others. To the extent that these 
single-inflaton estimations are reasonable, we again found plenty of roulette trajectories exist 
which satisfy the CMB+LSS data. Including the influence of the isocurvature degree of freedom 
on the observed power spectrum is needed to be completely concrete about such conclusions. 

Although there can be pure r-trajectories, we have not found pure 9 trajectories to be possible. 
The "natural inflation" model explored in [27j had a radial and an angular field, in a Mexican 
hat potential in which the Goldstone nature of the angular direction was broken into a cos(et2#) 
type term only after the radial motion settled into the (9-trough with random 9 values. (The 
a2 was also associated with non-perturbative terms similar to those invoked here.) The radial 
motion could come from either the small or large direction, the latter being like a chaotic inflation 
scenario for a first (unobserved) stage of inflation. The observed inflaton was to be identified with 
the angular motion from near a 9 maximum, where some fraction of the random 9 would reside, 
towards a 9 minimum. One of the features of the model was that near the maximum one could 
use the potential shape to get small gravity waves and yet significant scalar tilts. 

Why is our model so different? In a cos potential, inflation is possible in the immediate 
neighbourhood of the maximum, and so it is with our potentials if r could be fixed. However we 
have found that beginning in the 9 heights of our potentials near T^^min results in an outward 
radial flow, and, incidentally, insufficient e-folds to be of interest. These low r heights are not 
populated from inward r-flows. This is because angular dynamics are intimately tied to radial 
dynamics since the symmetry breaking of 9 and r are due to the same process, so hills and valleys 
in 9 are at all radii, ultimately guiding 9 trajectories to the r-valleys. And our approach to the 



25 



minimum is from a large radius flat potential rather than a growing chaotic inflation one. In 
general, we believe that our two-field Kahlcr modulus model gives a more natural inflation than 
natural inflation. 

The inflaton potential at large r is very shallow and in both the Kahler modulus and axion 
directions the motion rolls very slowly, even more so for larger r. In both r and 8 quantum kicks 
in the (quasi) de Sitter geometry can dominate over the classical slow roll drift. We identified 
this regime of the self-reproducing universe in our models. 

However, to have a self-reproduction regime depends on the immunity of the model against 
perturbative corrections, as we discussed in § [3] Corrections may lead to polynomial terms in 
t in V(r, 9) which may spoil the flatness, and possibly even spoil inflation itself. This problem 
is specific to the models of inflation associated with large t, as here. Models where inflation 
is realized with small values of the Kahler moduli, like "(better) racetrack inflation", are not 
sensitive to the alteration of the potential at large r. 

The issue is unclear because the exact form of the corrections is not yet known, so we are left 
with exploring the worst and the best case scenarios. The best case scenario is that corrections are 
absent or suppressed by a large volume factor I/V. The worst case scenario is where corrections 
generate significant terms for large r. We note that inflation based on the potential V(r, 8) with 
two variables is more protected from the corrections than inflation based solely upon the potential 
V(t). Indeed, if the potential V(r,8) is altered, the runaway character along the r-ridges may 
be changed to give a shallow minimum along r, and the region around this minimum might 
turn out to provide another suitable terrain for inflationary behaviour if there is slow roll in the 
8 direction. Although asymptotic flatness and self-reproduction may be destroyed by induced 
masses, the observable window at significantly smaller r may still allow inflations with enough 
c- foldings. As often happens in the investigation of string theory cosmology, assumptions such 
as the gentle nature of the uplift imposed here have to be made. That a viable stringy inflation 
seems feasible should motivate further work by the string theory community on the corrections 
and their role in inflation model building. 

In this paper, we did not treat the issues of reheating at the end of inflation and the relation 
of the model to the observed particle physics. We just assumed that all the inflaton energy is 
transferred to the energy of ultra relativistic particles of the Standard Model. We require that 
there should be no overproduction of dangerous particles such as other long-lived moduli or (non 
LSP) gravitinos. Such dangerous relics may "overdose the universe" or decay late enough to 
destroy the success of the Big Bang Nucleosynthesis of the light elements. Details of reheating 
depend on the Kahler moduli dynamics at the end of inflation and on the interactions of the 
Kahlcr moduli with other fields, including the Standard Model particles. The dynamics of t 
depends on the character of the potential around the minimum. In our picture all the other 
moduli except the inflaton Ti are stabilized and stay at the global minimum. If the value of 
T2,min exceeds the string scale, the QFT description of the inflaton potential around the minimum 
will be valid, with the inflaton beginning to oscillate after slow roll ends. The coupling to other 
degrees of freedom, in particular those of the Standard Model, that appear in the reheating 
process is unclear. The canonically-normalized inflaton may have interaction via the gravitational 
coupling and this interaction must be sufficiently suppressed (e.g., by the volume 1/V) to avoid 
significant radiative corrections to the inflaton potential. On the other hand, the oscillating 
inflaton should decay into Standard Model particles fast enough (in < 100 sec) to preserve 
successful Big Bang Nucleosynthesis. These requirements may put interesting constraints on the 
model for the couplings. 

An alternative possibility is that T2, m in is comparable to the string scale so that stringy effects 
play a direct role at the end of inflation. In this case one has to go beyond the QFT description of 
the processes. We note that small T2^ m in does not constrain the Kahler modulus/axion inflation, 
since the observed e-folds take place at large t: only the approach to, and consummation of, 
preheating would be affected, and would be quite different than in the QFT case. One may 
envisage the following as a possible scenario: the Kahler modulus, which corresponds to the 
geometrical size of the four-cycles, shrinks to a point corresponding to the disappearance of 
the hole. The energy of r cascades first into the excited closed string loops, then further into 
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KK modes in the bulk which interact with the Standard Model on the brane. Such a story is 
based on an analogy with the string theory reheating in warped brane inflation due to brane- 
antibrane annihilation investigated in [38 . In either the QFT or stringy case, reheating in Kahler 
moduli/axion inflation in the large volume stabilization model is an interesting and important 
question worthy of further study. 

Even within the context of the Kahler modulus/axion model explored here, the statistical 
element of the theory prior probability in the landscape of late stage moduli is unavoidable, so 
the terminology "roulette inflation" is quite appropriate — quite aside from the specific roulette 
trajectories we have identified that have a dominant angular motion before settling into a r- 
trough on the way to the minimum. Within the landscape, we may have to be content with error 
bars on inflationary histories that have a very large "cosmic variance" due to the broad range 
of theory models and trajectories possible as well as the data errors due to cosmic microwave 
background and large scale structure observational uncertainties. 
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